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Abstract 



We extract diffractive parton densities from data on diffractive deep inelastic 
scattering (DIS) and on diffractive photoproduction of jets. We explore the 
results of several ansatze for the functional form of the parton densities. Then 
we use the fitted parton densities to predict the diffractive production of jets 
and of VF's and Z's in pp collisions at the Tevatron. To fit the photopro- 
duction data requires a large gluon density in the Pomeron. The predictions 
for the Tevatron cross sections are substantially higher than data; this sig- 
nals a breakdown of hard-scattering factorization in diffractive hadron-hadron 
collisions. 
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I. INTRODUCTION 



In view of counterexamples [1,2] to the conjecture of factorization [3] of hard processes in 
diffractive scattering, it is important to test [4] factorization experimentally. In this paper, 
we present some results to this end. Specifically, we present fits^ to data from the ZEUS 
and HI collaborations on diffractive deep inelastic scattering (DIS) [6-8] and on diffractive 
photoproduction of jets [9]. Then we use these fits to predict cross sections in hard diffractive 
processes in pp collisions, with the assumption of factorization; we find that the predictions 
fail badly. 

We recall that diffractive events are characterized by a large rapidity gap, a region in 
rapidity where no particles are produced. We are concerned with the case where there 
is a hard scattering and where the gap occurs between the hard scattering and one of 
the beam remnants. Such hard diffractive events are observed in deep inelastic scattering 
(DIS) experiments [10] and are found to have a large rate: around 10% of the inclusive 
cross section. Diffractive jet production in pp collisions was earlier reported by the UA8 
collaboration [11], but under somewhat different kinematic conditions (larger There 
was also a report of diffractive bottom production [12]. Now, more diffractive data are 
being gathered from a variety of lepto-hadronic [6-9] and hadronic processes [13-18], but 
with substantially smaller fractions in the case of the diffractive production of jets and weak 
vector bosons in pp interactions than in DIS. 

Factorization for diffractive hard scattering is equivalent to the hard-scattering aspects of 
the Ingelman and Schlein model [3] , where diffractive scattering is attributed to the exchange 
of a Pomeron — a colorless object with vacuum quantum numbers. Ingelman and Schlein 
treat the Pomeron like a real particle, and so they consider that a diffractive electron-proton 
collision is due to an electron- Pomeron collision and that a diffractive proton-proton collision 
is due to a proton-Pomeron collision. Therefore they propose that diffractive hard cross 
sections are obtained as a product of a hard-scattering coefficient (or Wilson coefficient), a 
known Pomeron-proton coupling, and parton densities in the Pomeron. 

As was already known [2] before the advent of QCD, factorization is not expected to 
hold in general for diffractive hard processes. Furthermore, on the basis of a breakdown of 
the triple-Regge theory for soft single-diffractive excitation, Goulianos has proposed [19] to 
renormalize the Pomeron flux in an energy- dependent way. The agreement between data 
and his calculated cross sections is evidence that hard-scattering factorization is likely to 
break down in diffractive hadron-hadron collisions. 

However, one of us has recently proved factorization [20] for those diffractive hard pro- 
cesses that are lepton induced: these include diffractive DIS and diffractive direct photo- 
production of jets. The proof fails for hadron-induced processes. In this formulation, the 
primary non-perturbative quantities are diffractive parton densities [21-23] in the proton. 
Although we will use the terminology of "parton densities in the Pomeron", this mainly 
gives a useful way to describe the kind of parameterization we use for the diffractive par- 



^The fits presented in this paper represent a complete updating of our fits in an earlier preprint [5]. 

^By t we mean the invariant momentum-transfer-squared from the diffracted hadron. 
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ton densities, together with an indication of the quantum numbers that we behcvc to be 
exchanged across the rapidity gap. There is no necessary requirement that the object we 
call the Pomeron is the same as in soft scattering^. 

In principle, the parton densities in the Pomeron can be extracted from diffractive DIS 
(F2) measurements alone. Since the Pomeron is isosinglet and is its own charge conjugate, 
there is only a single light quark density to measure; one does not have the complications 
of separating the different flavors of quark that one has in the case of the parton densities 
of the proton. The Q dependence of the structure functions enables one to determine the 
gluon density. The HI collaboration has already presented [8] a fit of this kind. This type 
of data sufficiently determines the quark density in the Pomeron, and the HI fits suggest a 
large gluon content for the Pomeron. However, a more direct measurement of the gluons can 
be made in photoproduction, since the leading order processes have both quark- and gluon- 
induced terms. The ZEUS collaboration has already presented experimental evidence for 
a large gluon content of the Pomeron; they performed a combined analysis of their results 
on the diffractive structure function in deep inelastic scattering [6] and on diffractive jet 
photoproduction [9]. 

The main result of the ZEUS work was information on the overall normalization of the 
diffractive parton densities. In this paper, our aim is to obtain more detailed fits including 
the HI data, and to use the resulting fits to predict other cross sections. We use data on 
both DIS and photoproduction. Recently, the ZEUS collaboration has reported [25] new fits 
to their data that are made independently, but in a similar fashion to ours. 

For fitting the DIS data, we use full next-to-leading-order (NLO) calculations. The use 
of NLO rather than LO calculations is important since the gluon density is larger than the 
quark density. For the photoproduction data, we use leading-order calculations in a Monte- 
Carlo event generator in order to implement the experimental cuts. The event generator was 
constructed by two of us [26] as an extension to the POMPYT generator to allow the use 
of evolved parton densities in the Pomeron. With the resulting diffractive parton densities 
we cak:ulatc hard diffractive processes in hadron-hadron collisions, given the assumption of 
factorization. 

In the past, Ingelman and Schlein [3] and Bruni and Ingelman [27] have made similar 
calculations for one of the hadron-induced processes that we consider here {W/Z production). 
Their results have provided a commonly used benchmark in the phenomenology of these 
processes. They provide a choice of either "hard" or "soft" distributions of partons in the 
Pomeron, according to the /3 — > 1 behavior."^ The hard distributions give larger diffractive 
cross sections. At that time, there were no data to determine the distributions. We will find 
that although the quark distributions preferred by the DIS data are hard, our cross sections 
are substantially below those predicted by Bruni and Ingelman. We will present an analysis 
of the reasons for the lower values that we find. 

Nevertheless, our predictions for hadron-induced cross sections are well above the mea- 
surements [13-18], for both W production and jet production. In the case of jet production. 



^So Dokshitzer [24] would probably object to our use of the word "Pomeron" . 

■^Here, (3 is the fraction of the Pomeron's momentum that is carried by the struck parton. 
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the excess only occurs because of the large gluon density that is strongly preferred by the 
photoproduction data. 

This paper is organized as follows. In section II, we show our fits to diffractive deep 
inelastic and photoproduction data. In section III, we present some details of the formulae 
used to calculate the cross sections in hadron-hadron processes, and we discuss the kinemat- 
ics and phase-space cuts that we used. Then in sections IV and V, we present and discuss 
the results obtained for vector boson production and jet production, respectively. Finally, 
we summarize our findings in section VI. 

Other fits to the diffractive structure functions measured by HI have been made by 
Gehrmann and Stirhng [28] and by Kunszt and Stirling [21]. Golec-Biernat and Kwiecihski 
[29] assumed a parameterization of the parton densities in the Pomeron and found it to be 
compatible with the HI data on diffractive DIS. Their quark densities are about 30% smaller 
than ours, and they required the momentum sum rule to be valid. The new features of our 
work are a fit to a wider range of data, including photoproduction, the lack of an assumption 
of the momentum sum rule, and a calculation of the cross sections for diffractive jet and W 
and Z production, so as to test factorization by comparison with data from the CDF and 
DO experiments. 

II. PARTONS IN THE POMERON 

We will present a series of fits of parton densities in the Pomeron to data on diffractive 
DIS and diffractive photoproduction of jets. There are four sets of data that we use: 

• DIS data obtained by ZEUS using the rapidity gap method [6]; 

• DIS data obtained by ZEUS using their leading proton spectrometer (LPS) [7]; 

• DIS data obtained by HI using the rapidity gap method [8]; and 

• photoproduction data obtained by ZEUS using the rapidity gap method [9]. 



A. DIS 



Diffractive structure functions are related to the differential cross section for the process 
e + p ^ e+p + X: 



d/3dQ^dxpdt pQ^ 



where corrections due to exchange and due to radiative corrections have been ignored. 
Here xp is the fractional momentum loss of the diffracted proton (in the sense of light- 
cone momentum), and t is the invariant momentum transfer for the diffracted proton. The 
variables and y are the usual DIS variables, and /3 = x^j/x]?, with Xbj being the usual 
Bjorken scaling variable of DIS. 

Except for the ZEUS LPS data, the momentum transfer t is not measured, so we make 
fits to the structure function integrated over t, and write the structure function in the form: 
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F,''^'\p,Q^x^) = j\F?^^\p,Q\x^,t). (2) 

(We have set the lower hmit on t to —1 GeV^ to avoid including contributions where the 
putative diffracted proton results from fragmentation of a high pt jet. This point should 

not be important at small Xp. Moreover, the integrand in Eq. (2) is steeply falling in t so 
that the contributions to the integral from the region t < — 1 are quite small.) 

We next use hard-scattering factorization, proved in [20] , to write the diffractive structure 
function in terms of diffractive parton densities and hard-scattering coefficients: 

F^^^^ {p, Q^, Xp) = J2 ""ll^fa^^^ (/^' <3'' + NLO corrections, (3) 

a 

an equation valid to the leading power in Q. The hard-scattering coefficients are the same 
as in ordinary inclusive DIS. The predictive power of this equation comes from the DGLAP 
evolution equation obeyed by the parton densities and from the universality of the parton 
densities: they can be used to predict the cross sections for certain other diffractive hard 
processes. Factorization also holds for the diffractive structure function differential in t. 

We now assume that x^ is small. It is therefore sensible to use a parameterization of the 

dependence that is motivated by Regge theory. 

If Regge factorization is valid, then the dependence on x-p is of the form given by Regge 
theory, and therefore can be represented by a Pomeron flux factor, /p/p, that is related 
to the Pomeron-proton coupling mcasTircd in proton-proton elastic scattering. We do not 
necessarily expect Regge factorization to be valid. Nevertheless, we will assume that a 
suitable parameterization of the xp-dependence is of the Regge form, but possibly with 
different parameters than in proton-proton elastic scattering. If this form is not suitable, 
then we will find that we cannot fit the data, and a more general parameterization is needed. 
This can happen even though hard-scattering factorization remains valid in the form, (3), 
proved in Ref. [20]. 

So we will write the diffractive parton densities as a Pomeron flux factor times what are 
termed parton densities in the Pomeron: 

g^ xp) = /p/p(xp)/„/p(/?, (4) 

Furthermore we will assume the Pomeron flux factor is of the Donnachie-Landshoff (DL) [30] 
form: 

where is the proton mass, /5o — 1-8 GeV~^ is the Pomeron-quark coupling and ait) = 
Q;p-|-0.25i is the Pomeron trajectory. We treat ap as a parameter of our flts, instead of using 
the value given by Donnachie and Landshoff. Up to logarithmic corrections, the flux factor 
integrated over t is 



/p/p(xp) ^ cx;-^"^ (6) 

where C is a constant. 
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The t-dependence of the DL flux factor is not used in any of our fits, so the only use we 
make of the t dependence in Eq. (5) is to give a convention for a normahzation factor that 
is convenient for comparisons with other work. 

There is in fact another Pomeron flux factor that is commonly used, that of Ingelman 
and Schlein (IS) [3]. This differs from the DL flux factor primarily in its normalization. 
Since the same normalization factor appears in all our cross sections, its value is irrelevant 
to our phenomenology. Any change in the normalization factor is completely compensated 
by changing the parton densities by an inverse factor, and we obtain the parton densities 
from fltting a set of data without any a priori expectations as to their normalization. 

However, the normalization does affect the question of whether the momentum sum rule 
is obeyed by the parton densities in the Pomeron. Since it is not at present understood 
whether the sum rule is a theorem, this issue will not affect us. The momentum sum rule is 
not assumed in any of our flts. 

We will use Regge theory to make one further (correctable) assumption; this is, in effect, 
that the Pomeron has a deflnite charge conjugation parity and is an isosinglet. This implies 
that the distributions of u, d, u and d quarks are equal. Such an assumption is also valid in 
simple models where the rapidity gap is generated by gluon exchange. One possible mecha- 
nism for violation of the equality of the light parton densities would be the existence of an 
odderon, which has opposite charge conjugation to the Pomeron. The existence of Pomeron- 
odderon interference would break the equality of the quark and antiquark distributions. We 
will ignore this possibility, since there is no convincing phenomenological evidence to per- 
suade us of the odderon's existence. We also note that the issue does not concern us in 
DIS and photoproduction, since the hard-scattering coefficients are the same for quarks and 
antiquarks; in effect we will measure the average of the quark and antiquark distributions. 
Odderon contributions would only matter when we make predictions for diffractive cross 
sections at the Tevatron. 

In the data obtained using the rapidity gap method [6,8], the outgoing proton is not 
detected. Such data include "double-dissociative" contributions where the proton is excited 
to a state that escapes down the beam-pipe and thus, misses the detector. Factorization 
works for such final states, but since we will also wish to fit data where the outgoing proton 
is detected, we prefer to correct the data to remove the double-dissociative contribution. 
In the case of the ZEUS rapidity gap data [6], excited states up to about 4 GeV pass the 
diffraction selection cuts, and it is estimated that there is a contribution of (15 ± 10)% to 
the measured diffractive F2 from double-dissociative events. In obtaining our fits, we have 
corrected the relevant ZEUS data to take this into account. For the case of photoproduction 
data, we make the corresponding corrections for double dissociation and for nondiffractive 
contributions as well. No corrections have been made to the HI diffractive F2 data for 
which excited states up to 1.6 GeV are included. This point is relevant when we compare 
predictions obtained using our fits to data where the diffracted proton is detected (as in 
Sect. V) and also when we later compare our fits to both ZEUS and HI data. 



B. Photoproduction of Jets 

Similar formulae apply to photoproduction. For the direct diffractive photoproduction 
of a jet, 7 -|- p ^ jet -|- X -|- p, we let Et and r] be the transverse energy and pseudorapidity 
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of the jet. Then the cross section is 



= e/ dPh/,M /a/p(/^) ^"^g;;;" - (7) 



dx^dtdExdr] 



Here, da-y^a^j+x is the hard-scattering coefficient for the production of a jet in the colhsion 
of a photon and a parton of type a. It is the same as in inclusive photoproduction. The 
Pomeron flux factor /p/p and the parton densities in the Pomeron are the same as in the 
previous section. 

The proof [20] of the factorization theorem indicates that factorization is vahd for the 
direct part but probably not for the resolved part of diffractive photoproduction of jets. 
Fortunately, most of the cross section is from the direct process. This is known from the 
experimental data [9], and is also verified by our Monte-Carlo calculations. For the kinematic 
configurations of the data, we find a direct contribution that is 2 to 4.5 times larger than 
the resolved contribution, except at 77 = 0.75, where the two terms are comparable in size. 

If we use the factorization formula to calculate cross sections for diffractive photopro- 
duction, then presumably we should multiply the resolved contribution to the cross section 
by a correction factor^ similar to the one needed in hadron-hadron scattering (Sees. IV and 
V). Given the dominance of the direct contribution and the low precision of the current 
data — Fig. 3 — we feel that this is an unnecessary refinement in the present work. 

Beyond leading order, the separation between the resolved and direct processes is not 
unambiguous. Again, at the level of accuracy of the present data, we think that this is not 
an important enough issue to affect our results. 



C. Selection of data 

As far as the DIS data arc concerned, we restrict our attention to the subset of the data 
that is in the truly diffractive region. So we now explain the criteria we use to select the 
data for our fits. 

For the purposes of this paper, we define the diffractive component of a cross section to 
be the part of the cross section corresponding to the leading-power xp-dependence, of the 
form in Eq. (6). With this definition of diffraction, the cross sections reported by the ZEUS 
experiment [6, 7, 9] are the diffractive components. 

We do not need to address the question of whether the power dependence we use, with 
ttp around 1.1 or 1.2, is the ultimate asymptotic behavior as xp — > 0. We also do not 
require that this power be the same as in soft diffraction. It is sufficient that the power law 
represents an adequate approximation to a measurable part of the cross section, given that 
the factorization theorem [20] applies quite generally, and not just at small xp. This in fact 
implies that our restriction to diffractive data is mainly a matter of convenience, to reduce 



^On the basis of experimental evidence and of Regge models [2] for diffractive scattering, we might 
expect the correction factor to be less than unity, a suppression factor. However, the coherent 
Pomeron mechanism of CoUins, Frankfurt and, Strikman [1] would enhance the cross section. We 
will discuss this issue further in the conclusions. 
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the number of parameters, and to be in a region of xv common to the four sets of data to 
which we make fits. 

The HI data [8] include both non-diffractive and diffractive components, as is evidenced 
by the experiment's fit to their data with two powers of xp. To restrict our own fit to the 
diffractive region, we imposed the following cuts on all the DIS data: at /3 = 0.175 or 0.2, 
we require xw < 2 x 10~^, at (3 = 0.375 or 0.4, we require < 4 x 10~^, and a,t P — 0.65, 
we require ,xp < 1 x 10^^. We estimated these cuts by examining where the power-law 
associated with the Pomeron dominates Hi's fits to the xp dependence. The HI data at 
P — 0.1 and P — 0.04: are ehminated from our fits by this criterion. 

Another significant constraint is that we must restrict our fits to the truly deep-inelastic 
region. Outside of this region, the Icading-twist QCD factorization theorem for DIS does 
not hold. In fact some of the HI data lie very much in the resonance region. For example, 
they have points at /5 = 0.9 and = 4.5 GeV^. There the invariant mass of the excited 
hadronic system is mx = Q\J\I P — 1 = 0.7 GeV, i.e., close to the p resonance. While there 
are perturbative QCD results that apply in this region, they certainly do not include the 
usual inclusive factorization formula, Eq. (3). Most of Hi's data at /9 = 0.9 are at low mx-, 
while the data at smaller P have mx greater than about 2 GeV. Therefore we have simply 
chosen to discard the P = 0.9 data when we make our fits. 

With these cuts, the set of data which we fit comprises 77 points, of which 22 are from 
ZEUS DIS data obtained with the rapidity gap method [6], 3 from ZEUS DIS data from its 
LPS [7], 48 from HI DIS data [8], and 4 from the ZEUS photoproduction data [9]. These 
subsets of data we wiU call "ZEUS F2D3", "ZEUS LPS", "HI F2D3", and "ZEUS Photo", 
respectively. 

The region in which we make the fits is shown in Figs. 1, 2 and 3 below, which compare 
our fits to the data used to make the fits. 

D. Fits 

Each of our fits is represented by a parameterization of the initial distributions at 
Qq = 4 GeV^ for the u, u, d, and d quarks and for the gluon. The other quark distri- 
butions are assumed to be zero at this scale. For the DIS cross sections, we used full NLO 
calculations (with full evolution and with the number of fiavors set equal to 3), while for 
the photoproduction cross sections, we used a version of POMPYT that two of us have 
constructed [26], with the same evolved parton distributions as we used for DIS. The factor- 
ization and renormalization scheme is MS here and throughout this paper. As stated above, 
the Pomeron fiux factor is of the Donnachie-Landshoff form Eq. (5), but with an adjustable 
parameter for ccp, and we did not assume a momentum sum rule for the parton densities 
in the Pomeron, so that the choice of normalization for the fiux factor is irrelevant. The 
fits were made by minimizing x^, with the experimental systematic errors being added in 
quadrature to the statistical errors; no attempt was made to handle point-to-point correlated 
errors. The program used to perform the evolution was that of CTEQ [31]. 

We tried five functional forms for the parton densities, which we label A, B, C, D, and 
SG. For each of these five forms (Eqs. (8) and (9), below), we present the values of the 
parameters that give the best fit, given in turn each of the following three values of ctp: 
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ap = 1.08: 



Fit 


aq 


ag 


aq 


A 


0.496 ± 0.013 








B 


0.493 ±0.013 


9.3 ±2.5 





C 


0.501 ±0.022 





-0.008 ±0.031 


D 


0.565 ±0.026 


15.4 ±3.1 


-0.113 ±0.031 


SG 


0.470 ± 0.015 


12.6 ±3.9 





ftp = 


= 1.14: 






Fit 




ag 


aq 


A 


0.240 ±0.006 








B 


0.239 ±0.006 


4.5 ±0.5 





C 


0.249 ±0.011 





-0.031 ± 0.029 


D 


0.292 ±0.013 


9.7 ± 1.7 


-0.159 ±0.029 


SG 


0.225 ± 0.008 


7.4 ± 2.2 





ap = 


= 1.19: 






Fit 




ag 


aq 


A 


0.136 ± 0.004 








B 


0.135 ±0.004 


2.6 ±0.6 





C 


0.143 ±0.006 





-0.042 ± 0.028 


D 


0.175 ±0.008 


6.7 ± 1.0 


-0.191 ±0.026 


SG 


0.126 ±0.005 


5.0 ± 1.4 






TABLE I. Parameters of the fits for three different values of ap. 



• CKp = 1.08, which represents an appropriate value for a conventional Pomeron, as seen 
in soft scattering; 

• ap = 1.14, which approximates the best value of ctp associated with any of the param- 
eterizations except D; and 

• ar = 1.19, which gives the best fit associated with parameterization D. 

Since it is time-consuming to generate Monte-Carlo events for the photoproduction process 
and since the number of photoproduction data is small, wc first made some preliminary 
fits to DIS data alone to determine suitable values for ap, as listed above. Since the is 
not strongly dependent on ap, this seems to us to be sufficient. We will comment on the 
numerical values later. 

Four of the parameterizations, labeled A, B, C and D, use conventional shapes for the 
initial distributions. The final fit has a gluon distribution that is peaked near = 1, as 
suggested by the fit [8] obtained by the HI collaboration; we call this our "super-hard gluon" 
fit, SG. 
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Parameterizations A-D are all of the general form 

Pf,/AP, Ql) = s - /9) + 5, (1 - PT] , 

Pfg/AP,Ql) = agP{l-P), (8) 

with a series of constraints on the parameters. Note that since the Pomeron is isosinglet 
and sclf-charge-conjugate, the distributions of the u, d, u, and d quarks are all equal. Our 
first parameterization A, represents a conventional hard quark parameterization, where we 
set ciq = ttg — 0. Then in parameterization C we allow a soft quark term, while keeping no 
gluon term, so that Og = 0. In parameterization B we allow an initial gluon distribution but 
do not allow a soft quark term, so that dq — 0. Finally, in parameterization D we remove 
all the constraints. 

The super-hard gluon parameterization, SG, is of the form 

Pfq/AP,Ql)^aqP{l-P), 

Pfg/MQl)-a,p\l-(3f^ (9) 

i.e., the quark has a hard form, and the gluon is strongly emphasized at large p. The 
exponents for the gluon distribution were chosen somewhat arbitrarily, with no attempt 
being made to fit them. 

In Table 1, we show the parameters for each of the fits, and in Table 11 we show the 
values of x^, both for the total set of data and for each of the four subsets separately. In 
Figs. 1 and 2 we compare our fits to the HI and ZEUS DIS data. 

One important property of the fits is that the overall normalization of the quark distri- 
bution is quite well determined, as represented by the coefficient a^. This is not surprising, 
since the DIS cross section is dominated by a quark-induced process. The systematic shift 
to lower values of and Gg as ap increases is entirely due to the fact that the cross section 
has a factor l/xp"' and that the data are in the region < 10~^. 

The next important property is that a large initial gluon distribution is strongly preferred. 
This can easily be seen from the comparison to the photoproduction data in Fig. 3. With 
an initial gluon distribution that is zero, the cross section (dashed or dot-dashed curve) is 
over an order of magnitude below the data. Even though there are only 4 data points, the 
improvement when one goes to a parameterization with a large initial gluon distribution is 
the dominant effect in determining the gluon. The preference is also seen strongly in the 

values for the HI DIS data — see Table 11. However this preference is also associated 
with a negative soft-quark term in the initial parton densities, which would appear to be 
unphysical. We comment on this below. 

The relative size of the gluon distributions can be seen in Table HI, which gives the 
momentum sums for the initial parton distributions for each of the fits (in the case that ctp = 
1.14). Note that the total momentum sum, as opposed to its quark and gluon components, 
is invariant under DGLAP evolution. 

Now let us examine the fits in turn. 

Fit A has no gluons and no soft quark term. A good fit to the ZEUS DIS data is 
obtained: x^/(data point) is about 10/22 for the rapidity gap data and 2/3 to 5/3 for the 
LPS data (depending on the value of ap). However, only a moderately good fit is obtained 
for the HI data: x^/ (data point) ~ 70/48. The LPS data show a mild preference for a small 
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ap = 1.08: 



Fit 


Zs F2D3 


Zs LPS 


HI F2D3 


Zs Photo 


All Sets 


A 


8.2 


1.8 


81.9 


9.9 


101.8 


B 


5.9 


2.0 


77.7 


2.1 


87.8 


C 


8.5 


1.8 


81.6 


9.9 


101.8 


D 


9.3 


1.8 


65.3 


1.2 


77.6 


SG 


6.6 


1.9 


80.7 


2.1 


91.2 


«p = 


= 1.14: 










Fit 


Zs F2D3 


Zs LPS 


HI F2D3 


Zs Photo 


All Sets 


A 


8.6 


3.3 


68.8 


10.1 


90.8 


B 


5.8 


3.4 


65.3 


3.3 


77.8 


C 


9.7 


3.2 


66.8 


10.1 


89.7 


D 


10.8 


2.4 


42.1 


1.1 


56.3 


SG 


6.2 


3.7 


67.7 


1.9 


79.6 


ap = 


= 1.19: 










Fit 


Zs F2D3 


Zs LPS 


HI F2D3 


Zs Photo 


All Sets 


A 


9.5 


5.0 


72.3 


10.2 


97.1 


B 


6.3 


4.9 


68.3 


4.0 


83.6 


C 


11.0 


4.7 


69.0 


10.2 


95.0 


D 


12.4 


3.0 


34.8 


1.0 


51.2 


SG 


6.4 


5.8 


70.6 


2.0 


84.9 



TABLE II. for each of the fits. The data sets and the number of points are: ZEUS F2D3: 22 
points, ZEUS F2D3 LPS: 3 points, HI F2D3: 48 points, ZEUS photoproduction: 4 points. The total 
number of data points is 77. 



Fit 


Quarks 


Gluoii 


Total 


A 


0.160 





0.160 


B 


0.159 


0.750 


0.909 


C 


0.156 





0.156 


D 


0.133 


1.622 


1.755 


SG 


0.150 


0.375 


0.525 



TABLE III. Momentum sums 

/o d/5/?/a/p(/?) at Q = Qo = 2GcV for the fits with ap = 1.14. 
The quark column represents a sum over the 4 light flavors u, d, u, and d. 
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FIG. 1. Comparison of the fits for ap = 1.14 and the DIS data from HI [8] that were used in the 
fits. Fit A is represented by the dashed line, fit B by the dotted line, fit C by the dot-dashed line, fit D 
by the solid line, and fit SG by the heavy dashed line. 
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FIG. 2. Comparison of the fits for ap = 1.14 and the DIS data from ZEUS [6,7] that were used 
in the fits. The LPS data we used consist of just the three points at = 8GeV^. The code for the 
lines is the same as in Fig. 1. 
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FIG. 3. Comparison of the fits for ap = 1.14 and the ZEUS photoproduction data [9] used in the 
fits. The code for the lines is the same as in Fig. 1. 

value of ap, but this tendency is overwhelmed in the by a strong preference of the HI 
data for a larger value: ap = 1.14 gives a much better than ap = 1.08. However, the 
photoproduction data are not reproduced at all. 

Fit B differs from fit A by allowing an initial gluon distribution. Not surprisingly, this 
allows us to fit the photoproduction data much better with the x^/(data point) ranging 
from 2/4 to 4/4. However, these good values are mainly due to to the large errors in the 
last two data points. The dependence of the diffractive structure functions, as shown in 
Fig. 4, illustrates the strong influence of the gluon density on the evolution. We do not get 
a particularly good fit to the HI data. 

We next examine the effect of a soft quark term, in fit D. Although, in general, Regge 
theory suggests that there should be such a term in the parton densities at some level, 
what is surprising is that its coefficient is negative. The result is in fact the best of all 
our fits, including an excellent fit to the photoproduction data. The negative soft-quark 
term cannot represent the whole story, since it makes the initial quark densities negative 
at small j3. Notice that the quark distribution only becomes negative when f3 is below the 
region where we are fitting data, so that we do not have an unphysical quark density. If one 
wishes to extrapolate our parton densities to low /3 it would be sensible to replace the initial 
quark density by zero whenever the formula gives a negative value. This is in fact done 
automatically by the CTEQ evolution code that we are using, and one result of this can be 
seen in Fig. 5. In the curve for fit D at = 4.5 GeV^ there is a kink a little below (3 = 0.2. 
Notice that this kink disappears at higher Q^, when the effects of evolution give a larger 
positive contribution to the quark density at small (3. It is interesting that the restricted 
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FIG. 4. dependence of the fits with ap = 1.14. The code for the lines is the same as in Fig. 1. 
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FIG. 5. (3 dependence of the fits with ap = 1.14. The code for the lines is the same as in Fig. 1. 



set of data to which we fit provides no significant hint of a soft-quark term if we restrict to 
parton densities with no initial glue — as is seen in fit C. 

However, we are not sure to what extent the significance of this estimate of the soft-quark 
term is to be taken literaUy. If there were a systematic shift of the data with a series of 
points moving in a correlated way by about 1 standard deviation, the soft-quark term could 
be much reduced. Evidence that such a shift is possible is shown in Fig. 6. There we plot 
DIS data from both experiments. Generally the experiments are compatible, but there is a 
tendency for the ZEUS data to be about one standard deviation lower than the HI data for 
all the plots at /3 = 0.65. This would have a significant effect on the x^: at the level of ten 
units, given the number of points. 

Our final fits, SG, have an initial gluon density that is peaked at large /3, to mimic the 
one in the fits presented by HI [8]. Interestingly, we get a good fit to all but the HI data. It 
should be remembered, however, that we have found it appropriate to fit only to a subset 
of the data, as explained above, in Sec. II C. 

Finally, we comment on the value of ap. We find that we prefer the value 1.14 for 
Fits A, B, C, and SG. However, Fit D gives a value of 1.19. These values are certainly 
larger than the value for a soft Pomeron, and the lowest is given by fit D with ap = 
1.19, which is compatible with the value preferred by HI [8]: ctp = 1.203 ± 0.020 (stat.) ± 
0.013 (sys.) ^g°35 (model). However, observe that fit D with ctp = 1.14 provides a perfectly 
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FIG. 6. Comparison of the fits with the DIS data from both the HI experiment (closed circles) and 
the ZEUS experiment (open circles), in the region where both experiments have data. The code for 
the lines is the same as in Fig. 1. 
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adequate fit: = 56.3 for 73 degrees of freedom, and that the preference for the higher 
value of ap is entirely given by the HI data. 

In this context, it is worth examining Fig. 7, where we compare all the HI data with 
the predictions of our fits extended beyond the range where we make the fits. Some of the 
data are in a region of larger xp where we decided that the cross section is not Pomeron- 
dominated. The motivation for excluding certain regions of data can be seen particularly 
clearly in the graphs for j3 = 0.2. Furthermore at (3 = 0.9, the data appear to rise more 
steeply at small xp than our fits. While this is not conclusive, it suggests that a larger value 
of the Pomeron intercept, ap, would be needed to fit this subset of the data. As we explained 
in Sec. H C, these data are in the resonance region and it is thus not correct to include them 
in our fitting or to apply the factorization theorem in this region. 

Moreover, it has been established that the Pomeron trajectory is not universal, since the 
value of ap in hard scattering is not the same as in soft scattering. The proof of factorization 
certainly does not require such universahty. Therefore, there is no reason to assume that the 
same value of ctp applies to exclusive deep-inelastic processes and to the normal DIS region 
to which the factorization formula applies. 

Factorization does apply two constraints, however. First, parton densities are universal 
within the class of processes to which the factorization theorem applies, ap must be the 
same in these different processes. The second constraint arises from DGLAP evolution. 
Since evolution relates parton densities at different values of Q and /3, variations of ap with 
/3 and Q cannot be totally arbitrary. For example, suppose that at some particular value of 
Q, the value of ctp were larger at large f3 than at small Then evolution to larger Q would 
make the largest value of ap dominate at all /3. 

E. Shape of diffractive parton densities 

Since there are DIS data for a range of values of /3, the data do provide information on 
the shape of the diffractive quark distribution. For example, we are able to obtain significant 
information on the soft quark parameter in Eq. (8). 

However, we do not yet have similar information on the shape of the diffractive gluon 
distribution. We have two fits D and SG that provide good fits to the photoproduction data, 
but with dramatically different shapes. A direct measurement of the shape of the diffractive 
gluon distribution can be made in diffractive photoproduction of dijets by using the cross 
section da / djS'^^^ , where is the longitudinal (light-front) momentum of the jet pair 

relative to the Pomeron. In the leading-order parton approximation, j3'^^^ is exactly the 
momentum fraction of the parton in the Pomeron initiating the hard scattering. 

We see the implications of these observations in Fig. 8, where we superimpose our pre- 
dictions on preliminary data [25] for the diffractive photoproduction of dijets as a function 
of each of several kinematic variables. The only plot that enables us to distinguish the D 
and SG fits is that of the j3'^^^ dependence. The singular gluon is evidently preferred. In 
the other plots, both D and SG reproduce the normalization and shape of the cross section. 

In this same paper [25], fits of diffractive parton densities were presented. These were 
made independently of ours, with the same kinds of parameterization, but with the inclusion 
of only the ZEUS data; the results, particularly as regards the gluon distribution, are in 
general agreement with ours. 
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FIG. 7. Comparison of the fits for ap = 1.14 and all the DIS data from HI [8]. The code for the 
lines is the same as in Fig. 1. 
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FIG. 8. Comparison of the fits for ap = 1.14 and the preliminary 1994 ZEUS photoproduction 
data [25]. The code for the lines is the same as in Fig. 1. A double dissociation contribution of 
(31 zb 13)% has been subtracted from the data. Note that these data were not used in the fits. 

The preference for a singular gluon distribution is in agreement with the HI conclusions 
[8]. It is interesting that in our fits, the subset of the HI data that we use shows the opposite 
preference — see Table H. 



III. KINEMATICS AND CROSS SECTIONS FOR HADRON-INDUCED 

PROCESSES 

We now consider the production of W and Z bosons and of jets in diffractive pp collisions. 
In addition, wc consider W production with explicit measurement of the distribution of the 
final-state leptons. Schematically, these processes are 

p{Pi) + P{P2) (W or Z)+p + X, 
p{Pi) + P{P2) ^ jet + p + X 

p{pi)+p{p2)^{W ^l + v)+p + X. (10) 

We take the Pomeron to be emitted from the antiproton and the positive 2;-axis to be along 
the antiproton 's direction. 
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A. DifFractive jet production 



Consider the diffractive cross section for the production of a jet with rapidity y, in a 
hadron-hadron coUision. We will assume hard-scattering factorization [3,4]. At leading 
order (LO), the hard-scattering process is 2 —> 2 at the parton level, and results in a cross 
section of the form 

-^^Z2 J ^EtIEt J dy' j dXjpfw/p{xjp,ll) fa/p{Xa,l^) fb/F{Xb:l^)XaXb-^: (11) 
a,b 

where the sum is over all the active parton (quark, antiquark and gluon) flavors. The 
integration variables are E^, the transverse energy of the jet, y', the rapidity of the other 
jet, and Xp, the momentum fraction of the Pomeron. The momentum fractions of the 
partons, relative to their parent proton and Pomeron are 

x^^^(e-y + e-y') and Xb = ^^{e^ + e^'). (12) 

The functions fa/p{xa) and fb/f{xb) are the number densities of partons in the proton and 
Pomeron, respectively, while /p/p is the same Pomeron flux factor that we used in Sec. 11 A 
[see Eq. (5)]. da^^^/di is the partonic hard-scattering coefficient, and n is the factorization 
scale, which we set equal to Et- The hmits on the integrations are determined by the 
experimental conditions. 

The diffractive cross section given by Eq. (11) has the same structure as the factorized 
form of the corresponding inclusive cross section (i.e., without the diffractive requirement), 
except for the Pomeron ffux factor and the parton densities in the Pomeron. The same 
hard-scattering coefficient and nucleon parton distribution functions appear in both cross 
sections. 

The cross section given by Eq. (11) has contributions from a range of subprocesses. The 
indices a, b labeling the incoming partons range over the gluon and all the ffavors of quarks 
and antiquarks. The leading order (LO) form of the partonic cross section da-^^^ /di may be 
found in [32]. 

B. DifFractive W and Z production 

The cross section for the diffractive production of weak vector bosons is given by 

(7^^ = (Tq^ y2 [ — [ ^ [ — fp/p{xp) fb/v{Xb) fa/p{Xa) C^b^ ^ob ( , , (13) 

, J Xp J Xb J Xa \XaXbXf> / 

a,b 

where = \/2TTGFMy^/3s, Mys = Mw or Mz is the vector boson mass, Gp is the 
Fermi constant, Xh,Xa are momentum fractions of partons from the Pomeron and proton, 
respectively, and r = My^/s. For W bosons, C*^, = |Vgg/p if + e^/ = ±1 and zero 
otherwise, where q denotes a quark flavor, Cq the fractional charge of quark q and Vqqi is 
the Cabibbo-Kobayashi-Maskawa matrix element. For Z bosons, C*^ —1/2 — 2\eq\ siv? 9w + 
4|e5p sin^ Ow, where 6w is the Weinberg or weak- mixing angle. Similar expressions apply for 
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and Cgg which are relevant for gluon-induced scattering. The hard-scattering function 
cUab in the MS scheme and to next-to-leading order (NLO) in the QCD strong coupling ag 
can be found in [33]. 

C. DifFractive production of leptons from the W 

Since leptonic decays of W bosons include an unobserved neutrino, it is useful to compute 
the distribution of the observed charged lepton. The general formula for the distribution 
of leptons from W production has the same form as that for jet production, Eq. (11). In 
this case, we are only going to compute cross sections at leading order. Data have not been 
published for this particular cross section but since it is directly measurable, we think it is 
a useful quantity to work with. 

For the specific process — > {W~ — > e-|-Pe) + P + we have the leading-order cross 
section for quark-antiquark annihilation to a lepton pair: 



di 



(14) 



where Tw is the width of the W boson and u = —x^x^y/sETe ^. Using Eq. (14) in Eq. (11), 
one obtains the following cross section at the hadronic level: 



dy 



J 3;p 

a,b 



dEr fp/p{xv) fb/p{Xb) fa/p{Xa) 



QsTw[{Mw/2EtY - l]i/2 



(15) 



where Xa and are now given by 



Xa 



Xb = 



Mwe- 



M, 



w 



2Ei 



+ 



S XaXp 



'M, 



w 



2Ei 



- 1 



(16) 



We have suppressed the scale dependence of the functions Jijj in Eqs. (13) and (15); in 
actual computations, we set the scale equal to the vector boson mass. A similar equation 
may be obtained for the W'^ cross section. 



D. Inclusive cross sections 

Since we are particularly interested in the percentage of events that are diffractive, we 
also need to calculate the inclusive cross sections, that is, the ones without the diffractive 
requirement on the final state. The analog to Eq. (11) for the inclusive cross section for jet 
production is the standard formula 

(^^jet, incl ^ , , „ , , 



dy 

^ a,b 



/f da^ 
dET2ET j dy' fa/p{Xa, n) fb/piXb, IJ')XaXb—^, (17) 
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where Xa is given in Eq. (12), while Xb is now Xb = (e^ + e''^' ) Et / y/s ■ 

For the leptons from W~ production, the inclusive version of Eq. (15) is 



^^Icp, incl 

dy 



a.b 



w 

ab 



' 9^(9 



6srw[{Mw/2ETy - 1]V2 



(18) 



with a similar equation for W'^ production. In Eq. (18), u = —Xb\/sETe ^, defined 
in Eq. (16) while Xb is now given by Xb = M^/xaS. 

The analog to Eq. (13) for the inclusive total cross section for vector boson production 

is 



, J Xa J Xb 
a.b 



;VB , , 
ab ^ab 



. a, 



(19) 



IV. NUMERICAL CALCULATIONS OF W AND Z PRODUCTION 

For the calculations in this section, the factorization scale in the parton distributions was 
set to MvB- The values of the electroweak parameters that appear in the various formulae 
were taken from the particle data handbook [34], and we use only four flavors (u, d, s, c) in 
the weak mixing matrix, with the Cabibbo angle 9c = 0.2269. 



A. Comparison with previous calculations 

Bruni and Ingelman [27] computed diffr active W/Z cross sections neglecting any Q'^ 
evolution of the parton distributions in the Pomeron. At ^/s — 1800 GeV, they obtained 
the following diffractive fractions {R = a'^'^ /a''"'^): Rw++w- ^ 20% and Rz ~ 17% for total 
W and Z production, respectively. These rates are substantially larger than the few percent 
measured by CDF in [13]. 

As we will now explain, when one uses evolved Pomeron parton densities from our fits 
to data from HERA, one obtains substantially smaller rates than the Bruni-Ingelman ones. 
To understand these small rates, we first verify that we can reproduce the Bruni-Ingelman 
results. For these, we used their unevolved hard quark distribution in the Pomeron (given 
by their Eq. (4)), the same cut on xp: x^^"^ = 0.1, the EHLQl parton distributions in the 
proton and the Ingelman-Schlein (IS) flux factor:^ 

fw%M = J dt^ (3.19e«* + 0.2126^*) . (20) 

The flux in Eq. (20) differs by a factor of 1/2 from that in [27] because here we consider the 
case when only p diffracts while [27] considers the case when either p or p diffracts. 



^Note that since our purpose in using the IS flux is to compaxe our results with the Bruni-Ingelman 
calculations, we have used a Pomeron intercept of unity instead of the more accurate value used 
in our fits. 
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EHLQl 


EHLQl 


CTEQ4M 


Ref. [27] LO 


LO 


NLO 


LO 


NLO 


W+ + W- 


14000 


14300 


18100 


18700 


23500 


z 


1 100 


1100 


5600 


5500 


6900 



TABLE IV. Inclusive cross sections a^'^^^^^ (pb) for weak vector boson production. 

Next, wc evolved their Pomeron parton distributions and recalculated the cross sec- 
tions. Finally, to provide our best estimates of the rates, we repeated the calculations using 
CTEQ4M for the parton densities in the proton/antiproton and using our fits for the parton 
densities in the Pomeron. all with proper evolution. The cross sections were calculated using 
Eqs. (13) and (19) and the results obtained are summarized in Tables IV to VI. 

First, in Table IV, we show the inclusive cross section, o""^'^', which will give the denom- 
inator for the fraction of the cross section that is diffractive. We present the leading order 
(LO) result from [27] as well as our leading and next-to-leading order (NLO) results. At 
leading order, one observes that the use of the more up-to-date CTEQ4M densities in the 
proton instead of the EHLQl densities used by Bruni and Ingelman leads to cross sections 
that are 20% to 30% higher. Including the next-to-leading order contributions leads to 
another similar increase in the cross sections. 

The diffractive cross sections a^'^^^^ are shown in Tables V and VL In the columns 
labeled "BI" , we used the Bruni-Ingelman parton density in the Pomeron and the EHLQl 
parton densities in the proton, together with the Ingelman-Schlein flux factor (20). In the 
other columns we used our fits for the parton densities in the Pomeron together with the 
CTEQ4M parton distributions in the proton; we use the Donnachie-Landshoff form for the 
flux factor, (5), and ctp = 1.14. First, we use the same cut = 0.1 as was used by Bruni 
and Ingelman, to produce Table V. However, this allows to be rather larger than where 
Pomeron exchange is expected to dominate. So we also made calculations with = 0.01, 
for which the results are shown in Table VI. 

In column 3 of Table V we show our leading order results when we use the same unevolved 
parton densities as Bruni and Ingelman; we agree with their cross sections (column 2). Then 
we repeat the calculations but with correctly evolved parton densities in the Pomeron, with 
the Bruni-Ingelman formula being used as the initial data for the evolution at Qo = 4 GeV^ 
(column 5). The corresponding next-to-leading order cross sections are shown in columns 4 
and 6. We see that at either LO or NLO, evolution of the Pomeron parton densities leads 
to about a 30% reduction in the cross section. 

We also present in Table VII the diffractive fractions for total W production when either 
p or p diffracts. The fractions are obtained by dividing twice the diffractive cross sections in 
Tables V and VI (which are for single- sided diffraction) by the appropriate inclusive cross 
section in Table IV. 

The diffractive fraction obtained from the evolved BI Pomeron parton distribution, using 
columns 3 and 4 of Table IV for cr^.-^™^!^ jg about 14% for W production, compared with 
the 20% that is obtained using the unevolved BI Pomeron distributions. The corresponding 
percentages for Z production are a little smaller: 12% (evolved) and 17% (unevolved). 

In the last four columns of Tables V and VI we present the results when two of our 



24 



Pomeron: 


BI [27] 


BI 


BI 


Fit A 


Fit D 




unevolved 


unevolved 


evolved 


evolved 


evolved 


Proton: 


EHLQl 


EHLQl 


EHLQl 


CTEQ4M 


CTEQ4M 




LO 


LO 


NLO 


LO 


NLO 


LO 


NLO 


LO 


NLO 


W+ + W- 


1400 


1400 


1800 


1000 


1300 


300 


390 


650 


810 


z 


380 


380 


480 


260 


330 


77 


100 


170 


210 



TABLE V. DifFractive cross section cr^'^diff ^^^^ ^^j. ^ggi^ vector boson production when only 
p diffracts and with = 0.1. The cross sections using the BI distributions were computed with 
ap = 1, as in Ref. [27], but the cross sections using fits A and D were computed with ap = 1.14. 



Pomeron: 


BI 


Fit A 


Fit D 




evolved 


evolved 


evolved 


Proton: 


EHLQl 


CTEQ4M 


CTEQ4M 




LO 


NLO 


LO 


NLO 


LO 


NLO 


W+ + W- 


25 


38 


9 


14 


14 


21 


z 


3.2 


5.0 


1.1 


1.8 


1.6 


2.5 



TABLE VI. DifFractive cross section a^'^'^^f^ (pb) for weak vector boson production when only p 
diffracts, but now with x^^^ = 0.01. The cross sections using the BI distributions were computed with 
ap = 1, as in Ref. [27], but the cross sections using fits A and D were computed with ap = 1.14. 

fits (with ap = 1.14) shown in section II are used. Fit A is the one with a simple hard 
quark distribution and no glue as the initial values, while fit D, which has both quarks and 
gluons initially, is our best fit overall. Now fit A does not have the large gluon content 
that is necessary to fit the photoproduction data, so cross sections computed using fit A 
cannot be said to represent good predictions. However it is adjusted to fit DIS data, so that 
a comparison of predictions using fits A and D pinpoints situations where the large gluon 
distribution has a large effect. 

We have also calculated the cross sections resulting when we use the versions of the 
difFractive parton densities with a higher value of the Pomeron intercept, cip = 1.19. We 
find that the cross sections are reduced by 10% to 20%, depending on the value of a;™^. 
The reduced cross section arises because the diffractive parton densities are constrained to 
fit ZEUS and HI data at fairly small values of xp and we are now calculating cross sections 
at higher values of x^. So an increase in ctp results in a decrease in our calculated cross 
section for the hadron-induced processes. 

The LO and NLO cross sections resulting from fit A (columns 7 and 8 of Table V) 
are only about 30% of the evolved BI cross sections. (We will indicate the sources of this 
difference below, in Sec. IV B.) The diffractive fractions obtained from this fit, using the 
CTEQ4M entries in Table IV, are 3.3% (2.9%) for W (Z) production, as shown in Table 
VII. 

The quark distributions in fit D are about 20% higher than in fit A — see the values of 
Gq in Table I. However the cross sections for W production with fit D exceed those with 
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BI (unevolvcd) 
xg^ax ^ 0.1 


Fit A 
xg'^^ = 0.1 


Fit D 
xp^^^ = 0.1 


Fit A 
xg'^x ^ 0.01 


Fit D 
xg^ax ^ 0.01 


W+ + W- 


20% 


3.3% 


6.9% 


0.12% 


0.18% 


z 


17'X 


2.9'X 


(3.1'X 


0.05'X 


0.07'X 



TABLE VII. Diffractive fractions using NLO cross sections for + W and Z production when 
either p or p diffracts. 

fit A by a substantially larger factor, particularly at x^^^ = 0.1, where the cross section is 
more than a factor 2 higher. This arises because of evolution: the large gluon distribution 
in fit D increases the quark distribution at ^ = Mw compared with the case without the 
large gluon distribution. The increased quark density is most pronounced at small fractional 
momentum. Thus the effect is larger at Xp'^ = 0.1 than at Xp'^^ = 0.01, since in the first 
case, the quark from the Pomeron that makes the W has a smaller fractional momentum 
relative to the Pomeron. The NLO contributions further increase the fit D cross sections 
by 24%. Even so, the cross sections are still smaller, by a factor of 1.6, than the ones from 
evolved Bl Pomeron parton distributions. The rates from fit D (using NLO values) are 
6.9% (6.1%) for W (Z) production. These rates agree with those obtained by Kunszt and 
Stirling [21] with their Model 2 for diffractive parton distributions. 

The data from which our fits were extracted used a conservative cut on the Pomeron 
momentum, x'^^ — 0.01. The Pomeron fiux factor allows for the x-p dependence, but to 
ensure maximum compatibility with the HERA data without the assumption of standard 
Regge behavior, the same cut should be applied to the cross sections in hadron-hadron 
collisions. This results in the cross sections in Table VI, which therefore represent our 
most accurate prediction of diffractive W and Z production, given only the assumption of 
hard-scattering factorization, which of course we wish to test. Notice that with this cut 
the diffractive cross sections are over an order of magnitude smaller than with Xp'^^ = 0.1. 
The percentages obtained with this cut on xp for W (Z) production are 0.12% (0.05%) and 
0.18% (0.07%) using fits A and D, respectively, as shown in Table VII. The large reduction 
is due to the fact that we are not far from an effective kinematic limit: the cut on xp gives a 
maximum proton-Pomeron energy of 180 GeV, and partons typically do not carry the whole 
of the energy of their parent hadrons. 

B. Why are the fractions smaller than from BI? 

Although the data used in our fits support a "hard" quark distribution in the Pomeron, 
we predict that the diffractive W and Z cross sections are much smaller than those pre- 
dicted by Bruni and Ingelman, who also used hard quark distributions. For example, the 
diffractive fraction for W production computed using fit A, is six times smaller than Bruni 
and Ingelman's fraction (see Table Vll). 

Since Bruni and Ingelman's work served as an initial benchmark for subsequent work, 
it is interesting to understand the sources of this factor. We first address fit A, since that 
is our parameterization that is closest to Bruni and Ingelman's. The factor between the 
diffractive rates arises as an accumulation of several modest factors: 
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• A factor 0.9 because of the use of the CTEQ4M instead of the obsolete EHLQl distri- 
butions in the proton. (The denominator in the ratio of diffractive to inclusive cross 
sections increases by more than the numerator.) 

• A factor 0.7 for the effect of the evolution of the parton densities in the Pomeron. 

• A factor 1.7 for the use of the Donnachie-Landshoff flux factor instead of the Ingelman- 
Schlein flux factor, when the momentum sum is kept fixed. We have found that this 
factor arises from the following: 

- A factor of 2.5 to allow for our larger value of a^. 

— A factor of 0.7 to allow for the effects of the Pomeron slope a' . 

• A factor 0.16 because the DIS data indicate that the quarks have a momentum sum 
substantially less than the value of unity that was assumed by Bruni and Ingelman. ^ 

The first three factors in fact cancel. So one possible view is that the smallness of our 
results compared with those of Bruni and Ingelman arises essentially because of the change 
in the momentum sum of the quarks required by a fit to the data. An alternative view 
arises when one observes that the Pomeron-proton coupling is obtained by fitting data on 
high energy scattering, and that if one increases ctp, then the value of the Pomeron-proton 
coupling has to be decreased to keep the cross section at some particular energy fixed. (Of 
course, the energy dependence of the cross section would not be fitted so well.) The variable 
in high-energy scattering that corresponds to l/xf> for hard diffraction is s/M'^. Now, the 
typical value of xv in the data that we fit is about 10~^, which corresponds to s ~ 1000 GeV^, 
i.e., a fixed target energy of around 500 GeV. Therefore it is possible to argue that the factor 
of 2.5 for the larger value of ap should be combined with the factor of 0.16 for the momentum 
sum, to produce a factor of 0.4 for the momentum sum with the Pomeron-proton coupling 
fixed at a value appropriate for fixed-target CERN and Fermilab experiments. The overall 
reduction in our rates compared with those of Bruni and Ingelman then arises as a product 
of several factors, all less than unity. 

The effects of this decrease in cross section in going from the Bruni-Ingelman ansatz to 
our fitted distribution A are then somewhat compensated by the effects of the large gluon 
distribution we find in fit D. 

C. Lepton distributions for W production at the Tevatron 

In this section, we present our results for W production, but now with cuts on the emitted 
lepton /. Specifically, we calculate the electron's (or positron's) rapidity {y) distribution 
from Eq. (15) for the diffractive process, and Eq. (18) for the inclusive one. For the parton 
distributions in the Pomeron, we use our five fits with ap = 1.14 — see Eqs. (8)-(9) and 



^Note that in the case of the diffractive DIS cross section, this small momentum sum is mostly 
compensated by the effects of our larger value of ap, which increases the cross sections at small xp. 
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FIG. 9. Rapidity distribution of e~ in W~ production with the cut xp < 0.01. Solid curve results 
from using fit D, upper dashed with fit SG, dotted with fit B, lower dashed with fit A and dot-dashed 
with fit C. The lower dotted curve is the inclusive cross section scaled down by a factor 5 x 10"*^. 

Table I — evolved up to the W mass. We imposed a cut of 20 GeV on the Ex of the emitted 
lepton, and we integrated Xf up to — 0.01. 

Figure 9 shows our results for W~ production. For comparison, we also show the inclusive 
cross section rescaled by 5 x 10^'^ as represented by the lower dotted curve. The diffractive 
cross sections exhibit a strong fall-off in the region ye > —0.2 that is a consequence of the 
requirement of a rapidity gap. This fall-off is of course not present in the inclusive cross 
section. 

The diffractive cross sections are about 2% to 4% of the inclusive one at the left edge of 
the plots (at y = —3) depending on the fit used. At about y = —1.6 where the diffractive 
cross sections peak, this fraction drops to about 0.4% to 0.6% of the inclusive cross section. 
The cross sections using high glue fits D, B and SG, denoted by sohd, dotted and upper 
dashed curves, respectively, are larger (D gives largest cross section) than those using the low 
glue fits A (dashed) and C (dot-dashed) which overlap in the figure. The differences between 
the cross sections reffect first the size of the quark densities and then, in fits B, D and SG, 
the effects of a large gluon distribution on the evolution of the quark distribution. For 
example, fit SG has a smaller quark distribution than fit A, but its large gluon distribution 
pulls the cross section above that given by fit A. However the differences are moderate, at 
most a factor of 1.5. 

The corresponding cross sections for are shown in Fig. 10. The cross sections are 
larger than for the W^, because a valence up quark from the proton can be used to make a 
W~^, especially at large negative rapidities. In the plot, the rapidity gap exists for yg+ > —1.6. 



28 



dCT'"°'/dy* X (5E-4) 



FIG. 10. Rapidity distribution of e"*" in production with the cut xp < 0.01. The description of 
the various curves is the same as in the caption for Fig. 9. 

The same features as in the curves of Fig. 9 can be observed and thus, the same general 
inferences for W~ production can be made for this case as well. 



D. Comparison to CDF data for W production 

The CDF collaboration has presented data on diffractive W production from pp collisions 
at ^/s = 1800 GeV [13]. The W^s are produced with a rapidity gap in the region 2.4 < 
\r]\ < 4.2. They find that the fraction of diffractive to non- diffractive W production is [13] 
Rw = [1.15±0.51(stat) ± 0.20(syst)]%. This value corresponds to diffractive data corrected 
up to — 0.1 [35]. 

So, in Table VIII we present our diffractive fractions using Eq. (13) and our fits with 
ap = 1.14 for several different values of x'^^^. They are computed with the diffracted hadron 
being allowed to be either the proton or the antiproton. We see that for x™^^ = 0.01, the 
rates are an order of magnitude smaller while for x^^^ = 0.05, the rates are of the same 
order as the data. However, the preferred fits, with a large amount of initial glue (B,D and 
SG), yield rates which are about two to three times larger than the data. For x^^ = 0.1, 
our rates are a factor three to six larger than the central data value. 
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Fit 




x^ax = 0.05 


x^^^ = 0.1 


A 


0.12% 


1.9% 


3.3% 


B 


0.14% 


2.6% 


5.1% 


C 


0.12% 


1.8% 


3.2% 


D 


0.18% 


3.5% 


6.9% 


SG 


0.14% 


2.2% 


4.1% 



TABLE VIII. Diffractive fractions Rw for W production when either p or p diffracts. 

V. DIFFRACTIVE JETS 

In this section, wc present our results for jet production. We imposed the following cuts 
on the jet cross sections. These represent the effect of appropriate experimental cuts [14, 17] 
and of cuts to improve the significance of the signal. 

• We require that two jets are produced in the same half of the detector, i.e., yiy2 > 0, 
where Ui is the rapidity of jet i. This eliminates the region where the jets are in 
opposite hemispheres, since that region is well populated by non-diffractive events but 
is relatively unpopulated by diffractive events, because of the rapidity gap requirement. 

• Each jet is required to have a transverse energy Et greater than 20 GeV. This ensures 
that we are definitely in the perturbative region for the jets, but the cut could be 
relaxed. 

• Each jet's rapidity satisfies \y\ > Vcut = 1.8. 

Next, wc integrated over the rapidity of one of the jets to obtain a single jet distribution, 
but still subject to the above cuts on the other jet. Equations (11) and (17) were used 
for the diffractive and inclusive cross sections, respectively, with the parton distributions 
evolved to the scale Et- For the diffractive cross sections, the integral was performed up 
to = 0.01. In the following discussion, we will denote the rapidity of the final state jet 
by yjct instead of y. 

The resulting cross sections are shown in Fig. 11. There are no points in the middle part 
of the plot because of the rapidity cut. The cross sections using low glue fits A and C are 
nearly identical as depicted by the overlapping dashed (A) and dot-dashed (C) curves in 
the figure. The high glue fits D (solid curve), B (dotted) and SG (heavy dashed) yield cross 
sections that are about an order of magnitude larger, with D being largest, than those using 
low glue fits. This difference reflects the sensitivity of this particular type of cross section 
to the gluon content of the Pomeron. The lower dotted curve, which is symmetric about 
y = 0, represents the inclusive cross section scaled down by a factor of 5 x 10~^. 

The diffractive jet percentages are shown in Fig. 12, where i?x 100 is plotted as a function 
of yjet, with R = ^^jct.inci/|j^''.°|^ ■ One finds that the rates R are largest when fit D is used, 
varying from 2.7% to 5.7%. With fit B, whose gluon distribution is about a factor of two 
lower than fit D, the rates are also about a factor of two smaller. The rates obtained with 
fits A and C are much lower ranging from 0.2% to about 0.3%. With fit SG, the resulting 
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FIG. 11. Rapidity distribution of jet cross sections, with Et > 20GeV, y > 1.8, and Xf < 0.01. 
The description of the various curves is the same as in the caption for Fig. 9. 

curve is relatively flat giving a rate of about 1.2%. The rates are largest at j/jet = —4, 
then decrease as yjet increases. Of course, the large rates for distributions D, B and SG, all 
with the large gluon distribution, directly result from the fact that there is a gluon-induced 

subprocess. 

We end this section by making comparisons with data on diffractive dijet production 
from CDF and DO at y/s = 1800 GeV. CDF has measured dijet data both with a rapidity 
gap requirement [14] and with Roman pots [15] along the antiproton beam direction. In 
the first case, the cross section for dijets produced opposite a rapidity (r)) gap in the region 
2.4 < \ri\ < 4.2 is measured. Each jet is required to have a minimum of 20 GeV and 
rapidity 1.8 < \ri\ < 3.5. They also measure the dijet cross section without a rapidity gap, 
i.e., what wc refer to in this paper as the inclusive cross section. The diffractive fraction they 
measure is [14] Rjj = [0.75 ± 0.05(stat) ± 0.09(syst)]%. This measured value is appropriate 
for Xf> < 0.1 [35]. The fractions that we obtain using the above cuts and our fits with 
ctp = 1.14 are shown in Table IX, for several values of x^'^^. Our calculation assumes that 
either the antiproton or the proton is diffracted. The rates obtained with fit D or B are 
from 3 to 22 times larger, while those obtained with fit C or A range from being about 70% 
smaller to being a few percent larger than the measured value, depending on the value of 
Xp^. The rates using fit SG are also significantly greater than the data but smaller than 
the rates with fits B and D. This reflects the low number of gluons in flt SG; they were more 
effective in the photoproduction at producing jets because of their relatively large fractional 
momentum relative to the Pomeron. 
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FIG. 12. Diffractive jet production percentages, with Et > 20GeV, y > 1.8, and xp < 0.01. Solid 
curve results from using fit D, upper dashed with fit SG, dotted with fit B, lower dashed with fit A and 
dot-dashed with fit C. 



Fit 






x'p^ = 0.1 


A 


0.23% 


0.63% 


0.83% 


B 


1.9% 


6.0% 


8.1% 


C 


0.23% 


0.61% 


0.79% 


D 


3.9% 


12.3% 


16.4% 


SG 


1.2% 


2.5% 


3.2% 



TABLE IX. Diffractive fractions Rjj for dijet production when either p or p diffracts and using 
cuts on Et and y appropriate for the CDF rapidity gap data. 
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Fit 


A 


B 


C 


D 


SG 


Rjj 


0.20% 


1.8% 


0.19% 


3.7% 


0.85% 



TABLE X. Diffractive fractions Rjj for dijet production when only p diffracts and using cuts on 
XV, Et and y appropriate to the CDF Roman pot data. 



Fit 


A 


B 


C 


D 


SG 


Rjj 


0.59% 


5.8% 


0.57% 


11.8% 


2.4% 



TABLE XL Diffractive fractions Rjj for dijet production when either p or p diffracts and using 
cuts on xp, Et and y appropriate for the DO data. 

With their Roman-pot-triggered diffractive sample, CDF has measured a diffractive frac- 
tion of Rjj = [0.109 ± 0.003 ± 0.016]%. The data in this sample correspond to xp in the 
range 0.05 < Xp < 0.1, with the jets having minimum of 10 GeV. The fractions we obtain 
using the same kinematic cuts and our fits with ap = 1.14 are presented in Table X. In this 
case, our calculation assumes that only the antiproton is diffracted. The ones obtained with 
fits D, B and SG are from 8 to 34 times larger than the data, while those obtained with fits 
C and A are about twice as large. 

Finally, DO also has some preliminary data [17] on diffractive dijet production. They 
require a rapidity gap opposite the dijets, which have = 12 GeV and |?7jet| > 1-6. The 
diffractive fraction they measure with an estimated = 0.03 is Rjj — [0.67 ± 0.05]%. 
Our calculated fractions are shown in Table XI; as with our previous calculations, we use 
the fits with ap = 1.14 and assume that either the antiproton or the proton is diffracted. 
The realistic fits (with a large gluon content) are well above the data, by factors of 9, 18, 
and 4 for fits B, D, and SG, respectively. The cross sections obtained from fits A and C 
are a bit smaller than the data; these fits give a correct normalization for diffractive DIS, so 
again we see the importance of the photoproduction data in demonstrating a breakdown of 
factorization. 

VI. CONCLUSIONS 

We have presented parton distributions in the Pomeron resulting from fits to data on 
diffractive DIS and diffractive photoproduction at HERA. In order to explore the sensitivity 
of the data to different aspects of the parton densities, we made several fits with different 
assumptions for the initial parton densities. We find that only those parameterizations with 
a large amount of ghic (B, D and SG) are able to provide a good fit to the photoproduction 
data. The other two parameterizations (A and C), which are constrained to have no gluons 
in the starting distributions, badly underestimate the photoproduction cross sections. 

We also find that the normalizations of both the quark and gluon densities in the Pomeron 
are well determined by the data. As regards the shape, hard distributions are preferred. 
But in the case of the gluon, the question still remains as to whether a conventional hard 
distribution (1 — /3 at large /3) or something harder is correct. We are able to obtain 
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satisfactory fits with both a hard gluon, in fit D, and a harder gluon, in fit SG. We have 
shown how measurements of the f3 dependence of the photoproduction cross section will be 
able to provide much better information. 

Prom our fits, we predicted the cross sections for vector boson production and dijet 
production in diffractive pp interactions at the Tevatron. The rates represent a realistic 
prediction of the cross sections, given the assumption of factorization. We find that the 
predictions are a factor of several above the measured cross sections. In the case of the 
jet cross sections, it is only for the physically correct "high-gluon" fits that the predictions 
substantially exceed the data. The lack of agreement between the predictions and the data 
indicates a substantial breakdown of factorization in diffractive pp interactions. 

For the predictions to match the measured diffractive rates of W production by CDF, 
suppression factors (= prediction/data) ranging from three to six must be applied. In the 
case of diffractive dijet production, the suppression factors appear to be somewhat larger, 
around 10. (We refer only to the realistic fits, with a large amount of glue.) 

Further work to measure the suppression factors is necessary to obtain a fuller under- 
standing of the dynamics of diffractive hadron-hadron interactions. One interesting possi- 
bility is to search for the contribution predicted by the coherent Pomeron mechanism of 
Collins, Frankfurt and Strikman [1] , which in fact gives an enhancement of the cross section 
at large Such an enhancement is suggested by the UA8 data [11]. This and our results 
on photoproduction show that the measurement of /3 distributions is important. It should 
be noted that the UA8 data are at larger \t\ than the data which we have fitted. 

With regards to extracting diffractive parton densities, further work is also needed to 
understand the differences between the ZEUS and HI data, as illustrated in Fig. 6. The 
differences arc suggestive of a systematic error that is correlated point-to-point. This indi- 
cates that we need to be careful about taking the values at face value, and in fact that 
systematic errors need to be treated more correctly. Our negative value for the soft quark 
term in fit D is worrying; note that it is driven by the HI data, as can be seen from the 
vahics in Tabic 11. 

It is also important to test the universality of ap, for example, to test whether its value 
is different in exclusive and inclusive processes, as is suggested by Fig. 7. 

Finally, further tests of factorization can be accomplished at HERA. For example, we 
expect hard-scattering factorization to be valid for heavy quark production in DIS as well, 
but not for any resolved photoproduction process. 

Note added: After completion of the work for this paper, a paper by the ZEUS col- 
laboration [36] appeared. It provides the official version of the diffractive photoproduction 
data [25] that we showed in Fig. 8, and the paper reports a QCD analysis of the ZEUS data. 
This analysis was performed independently of the one in the present paper, but in a similar 
style, and the conclusions as regards the parton densities in the Pomeron are similar. 
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